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ABSTRACT 

The Stardust mission was successfully launched on 
February 7. 1999. It will be the first mission to re- 
turn samples from a comet. The sample return cap- 
sule, which is passively controlled during the fastest 
Earth entry ever, will land by parachute in Utah. The 
present study describes the analysis of the entry, 
descent, and landing of the returning sample capsule 
utilizing the final, launch configuration capsule mass 
properties. The effects of two aerodynamic insta- 
bilities are revealed (one in the high altitude free 
molecular regime and the other in the tran- 
sonic/subsonic flow regime). These instabilities 
could lead to unacceptably large excursions in the 
angle-of-attack near peak heating and main para- 
chute deployment, respectively. To reduce the ex- 
cursions resulting from the high altitude instability, 
the entry spin rate of the capsule is increased. To 
stabilize the excursions from the transonic/subsonic 
instability, a drogue chute with deployment trig- 
gered by a gravity-switch and timer is added prior to 
main parachute deployment. A Monte Carlo dis- 
persion analysis of the modified entry (from which 
the impact of off-nominal conditions during the en- 
try is ascertained) predicts that the capsule attitude 
excursions near peak heating and drogue chute de- 
ployment are within Stardust mission limits. Addi- 
tionally, the size of the resulting 3-a landing ellipse 
is 60.8 km in downrange by 19.9 km in crossrange, 
which is within the Utah Test and Training Range 
boundaries. 

NOMENCLATURE 

CFD Computational Fluid Dynamics 

c.g. Center-of-Gravity 

DOF Degree-of-Freedom 

DSMC Direct Simulation Monte Carlo 


LAURA 

Langley Aerothermodynamic Upwind 
Relaxation Algorithm 

NASA 

National Aeronautics and Space 
Administration 

NNE 

North-Northeast 

PICA 

Phenolic Impregnated Carbon Ablator 

POST 

Program to Optimize Simulated 
Trajectories 

SRC 

Sample Return Capsule 

ssw 

South- Southwest 

TLNS3D 

Thin Layer Navier Stokes 
Three-Dimensional 

UTTR 

Utah Test and Training Range 

a T 

Total angle-of-attack (angle between the 
spin axis and atmospheric velocity 
vector), deg 

Y 

Flight-path angle, deg 


INTRODUCTION 


On February 7, 1999, the fourth of NASA’s 
Discovery class missions, known as Stardust, was 
successfully launched from Cape Canaveral Air Sta- 
tion. Stardust is a comet sample return mission. It 
will be the first mission to return samples from a 
comet. The spacecraft will encounter the comet 
Wild-2 in 2004. Stardust will come within 100 km 
of the comet nucleus and deploy a sample tray to 
collect cometary and interstellar dust particles 
(Fig. 1). Upon Earth return in January 2006, the en- 
try capsule (Fig. 2), containing the comet samples, 
will be released from the spacecraft and land by 
parachute at the Utah Test and Training Range 
(UTTR). The entry velocity will be the highest of 
any Earth-returning mission (inertial velocity about 
12.9 km/s). A new heat shield made of PICA (Phe- 
nolic Impregnated Carbon Ablator) will be used to 
protect the Sample Return Capsule (SRC) from the 
intense heat of reentry. 1 
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Figure 1. Stardust spacecraft flight configuration. 



Figure 2. Stardust sample return capsule 
configuration. 


Four hours prior to entry, the SRC will be spun- 
up and separated from the main bus. The SRC has 
no active control system, so the spin-up is required 
to maintain its entry attitude (nominal 0 deg angle- 
of-attack) during coast. Throughout the atmospheric 
entry, the passive SRC will rely solely on aerody- 
namic stability for performing a controlled descent 
through all aerodynamic flight regimes: hypersonic- 
rarefied, hypersonic-transitional, hypersonic- 
continuum, supersonic, transonic, and subsonic. The 
SRC must possess sufficient aerodynamic stability 
to overcome the gyroscopic (spin) stiffness in order 
to minimize any angle-of-attack excursions during 
the severe heating environment. Additionally, this 
stability must persist through the transonic and sub- 
sonic regimes to maintain a controlled attitude at 
parachute deployment. 

This study describes the analysis of the entry, 
descent, and landing sequence for the returning 
sample capsule. This analysis consists of performing 
a trajectory simulation of the entire entry (from bus 
separation to landing) to predict the descent attitude 
and landing conditions. In addition, a Monte Carlo 
dispersion analysis is performed to ascertain the 
impact of off-nominal conditions which may arise 


during the entry to determine the robustness of the 
Stardust SRC design. Specifically, the SRC attitude 
near peak heating and parachute deployment is of 
interest, along with the landing footprint ellipse. 

The SRC is restricted to land within the UTTR 
site. For mission success, a high-fidelity aerody- 
namic database is essential for accurately predicting 
the landing location, as well as the attitude of the 
SRC at critical phases (e.g., peak heating and para- 
chute deployment) during the entry. In this paper, 
the aerodynamics utilized in the entry simulation is 
discussed first, followed by a description of the 
nominal entry sequence of the SRC. Finally, the 
results of the Monte Carlo entry dispersion analysis 
are presented. 

ANALYSIS 

Aerodynamics 

The aerodynamic database utilized for the SRC 
in the flight simulation studies is constructed from a 
combination of computational fluid dynamics (CFD) 
calculations, and wind tunnel and existing historical 
flight data as described by Mitcheltree, et. al. 2 This 
large variety of sources for the aerodynamics is re- 
quired because the SRC traverses many different 
flow regimes (hypersonic-rarefied, hypersonic- 
transitional, hypersonic-continuum, supersonic, tran- 
sonic, and subsonic) during its entry. At the outer 
reaches of the atmosphere, free molecular flow 
aerodynamics are employed. In the rarefied flow 
regime, Direct Simulation Monte Carlo (DSMC) 
simulations are used to define bridging functions for 
the aerodynamic coefficients. In the hypersonic- 
continuum regime, a matrix of solutions from the 
computational fluid dynamics code LAURA (Lan- 
gley Aerothermodynamic Upwind Relaxation Algo- 
rithm) 3 describe the aerodynamics. At supersonic 
and transonic speeds, the aerodynamics are taken 
from two sets of existing wind tunnel data, aug- 
mented by CFD results from the TLNS3D (Thin 
Layer Navier Stokes Three-Dimensional) code 4 , 
along with ballistic range tests. Subsonic aerody- 
namics are defined by a combination of static wind 
tunnel measurements and dynamic free flight meas- 
urements. 5 These sources are blended to form a 
comprehensive database which describes the aero- 
dynamics of the SRC for the expected flight condi- 
tions. Figure 3 shows the range of application of the 
various aerodynamic sources mentioned above. The 
aerodynamic characteristics of the SRC are de- 
scribed in detail in Ref. 2. 
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Figure 3. Stardust SRC aerodynamic database. 

Trajectory Simulation 

The trajectory analysis is performed using the 
six- and three-DOF (degree-of-freedom) versions of 
the Program to Optimize Simulated Trajectories 
(POST). 6 This program has been utilized previously 
for similar applications. 7-9 The three-DOF program 
(which integrates the translation equations of mo- 
tion) is used from bus separation to atmospheric 
interface. The six-DOF version of POST (which 
integrates the translational and rotational equations 
of motion) is used from atmospheric interface until 
parachute deployment. The three-DOF program is 
used again from parachute deployment to landing. 
The trajectory simulation includes Earth atmos- 
pheric (GRAM-95) 10 and gravitational models, cap- 
sule separation and non-instantaneous parachute 
deployment models, and capsule aerodynamics and 
mass properties. The validity of the present approach 
has been demonstrated recently through compari- 
sons between the Mars Pathfinder pre-flight predic- 
tions of the flight dynamics and the flight data. 11 

During the entry, off-nominal conditions may 
arise which affect the descent profile. These off- 
nominal conditions can originate from numerous 
sources, such as capsule mass property measurement 
uncertainties, separation attitude and attitude rate 
uncertainties, and limited knowledge of the flight- 
day atmospheric properties (density, pressure, and 
winds). Additionally, computational uncertainty 
with the aerodynamic analysis and uncertainties with 
parachute deployment are contributing sources of 
uncertainty. In this analysis, an attempt is made to 
conservatively quantify and model the degree of 
uncertainty in each mission parameter. For this mis- 
sion, 41 potential uncertainties were identified. 
These uncertainties are grouped into two categories 
(exo-atmospheric and atmospheric) and are listed in 
Tables 1 and 2, respectively, along with the corre- 


sponding 3-a variances. For most of the parameters, 
a Gaussian distribution is sampled. Flowever, for the 
center-of-gravity (c.g.) offset quadrant and parachute 
deployment parameters (gravity-switch, timers, and 
aerodynamics), uniform distributions are utilized to 
model their operating performance. 


Table 1. Exo-Atmospheric Mission Uncertainties 


Mass Properties 

3-a Variance 

Mass 

. +4.7e-5 kg 

eg position along spin axis. 

.±0.0254 cm (0.01 in.) 

eg position off spin axis 

. ±0.0254 cm (0.01 in.) 

Major moment of inertia 


C-^-xx’ Iyy> ^zz) 

. ±2%, ±5%, ±5% 

Cross products of inertia 


(Ixy’ -^xz’ lyz) 

. ±0.003 kg-m 2 

Post-Separation State Vector 

Position 1 correlated with coveriance 

Velocity/ matrix corresponding to Ay : = ±0.055 deg 

Pitch/yaw attitude 

. ±0.56 deg 

*Pitch/yaw rate 

. ±2.828 deg/s 

Roll rate 

. ±2 rpm 

Separation 

Spring induced AV : 


Body x-axis AV 

. ±0.0508 m/s (2 in/s) 

Body y-axis AV 

. ±0.0097 m/s (0.38 in/s) 

Body z-axis AV 

. ±0.0097 m/s (0.38 in/s) 


*Uncertainty sampled using uniform distribution 


As will be shown in the results, the successful 
return of the cometary samples by the Stardust SRC 
depends heavily on the validity of the Monte Carlo 
analysis. Increased reliance on entry simulations for 
mission success places considerable importance on 
selecting appropriate uncertainties. As confidence 
increases in the analysis accuracy, cheaper and/or 
higher performance entry systems can be selected 
for future missions. 

RESULTS AND DISCUSSION 
Nominal Mission 

Original Entry Sequence 

In the original nominal Stardust entry sequence, 
the SRC enters the atmosphere with a spin rate of 5 
revolutions per minute (rpm). The spin rate main- 
tains entry attitude (nominal 0 deg angle-of-attack) 
until atmospheric interface (since the SRC possesses 
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Table 2. Atmospheric Mission Uncertainties 


Aerodynamic 

3-a Variance 

Free molecular aerodynamics, 


c A 

.±10% 

Cv c Y 

.±8% 

C m , C 

.±12% 

Hypersonic continuum aerodynamics. 

c A 

.±4% 

C N , C Y 

.±8% 

C m , C 

.±10% 

Supersonic continuum aerodynamics, 

c A 

.±10% 

Cv c Y 

.±5% 

C m , C 

.±8% 

Subsonic continuum aerodynamics, 


c A 

.±5% 

*Hypersonic dynamic stability coefficients. 

C C 

> “'mq’ v -'nr 

.±0.28 

♦Supersonic dynamic stability coefficients. 

C C 

v -'nr 

.±0.2 

Atmosphere 

Pressure, density, winds: 


GRAM-95 model 

. 3-a scale factor 

Other 

Ablation mass 

.±10% 

♦Drogue chute g-switch 

.±10% 

♦Droaue chute deployment timer 

.±0.05 sec 

♦Drogue chute aerodynamics, 

.±10% 

♦Main chute deployment timer 

.±0.05 sec 

♦Main chute aerodynamics, C ; \ 

.±10% 


*Uncertainty sampled using uniform distribution 


no active control system). As the SRC descends, it 
must rely solely on atmospheric stability in all flow 
regimes to minimize any angle-of-attack excursions 
until main parachute deployment at Mach 0.16. 
However, an analysis of the entry capsule revealed 
that the SRC was statically unstable in the free mo- 
lecular flow regime due to its aft center-of-gravity 
location (0.283 m or 0.349 body diameters back 
from the nose). 2 The six-DOF analysis reveals that 
the pitch rate induced by the instability during the 
free molecular regime carries into the transitional 
region where high angles-of-attack are produced. 
This static instability causes the SRC to pitch-up to a 
total angle-of-attack (a T ) above 70 deg as seen in 
Fig. 4. As the SRC descends into the continuum 
regime (where it is statically stable), the angle-of- 
attack damps out and reduces to small values. 



Time, sec 


Figure 4. Original nominal Stardust attitude profile. 


High angles-of-attack early in the entry are a 
concern, since they can lead to angles-of-attack 
greater than 10 deg at peak heating (which is a Star- 
dust mission limit). Angles-of-attack greater than 10 
deg increase afterbody heating near the shoulder 
regions and can damage the afterbody thermal pro- 
tection system. Moreover, high angles-of-attack are 
worrisome because the SRC is stable flying back- 
wards. Off-nominal attitude and attitude rate condi- 
tions at atmospheric interface could result in a 
backward entry attitude leading to a loss of the cap- 
sule. In addition to the high altitude static instability, 
a low altitude dynamic instability in the transonic 
and subsonic flow regimes also exists, again caused 
by the aft center-of-gravity location. This tran- 
sonic/subsonic dynamic instability could induce a 
tumbling motion prior to main parachute deploy- 
ment, raising concerns for a successful deployment. 
Reference 2 describes the SRC's stability/instability 
in the various aerodynamic regimes in more detail. 

Modified Entry Sequence 

The original nominal entry sequence is unac- 
ceptable since off-nominal conditions could result in 
mission failure. Modification of either the entry se- 
quence or the SRC is required to improve the prob- 
ability of mission success. 

The high altitude and transonic/subsonic insta- 
bilities could be eliminated by moving the c.g. of the 
SRC forward to 0.26 body diameters back from the 
nose. 2 Since the SRC has only three major compo- 
nents (forebody heat shield, sample tray, and para- 
chute canister), movement of the c.g. via re- 
packaging of the SRC is difficult. The size and mass 
of the sample tray preclude large movements in the 
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c.g. Ballast could be added to the nose of the SRC to 
move the c.g. forward. However, a prohibitive large 
amount (22.4 kg) of ballast is required to move the 
c.g. sufficiently forward to remove the instabilities. 

Augmentation of the SRC’s stability is required 
to eliminate the large angle-of-attack excursions. 
Several concepts were considered. For example, 
adding an aft-skirt would provide a restoring torque 
in the free-molecular regime stabilizing the SRC. 12 
However, these devices, once serving their purpose, 
must be discarded in order to avoid destabilizing the 
lower flight regimes. To avoid such complications, 
the solution selected to address the high altitude in- 
stability is to increase the spin rate of the SRC upon 
entry. The higher spin rate, although not eliminating 
the instability, increases the gyroscopic stiffness of 
the SRC sufficiently to retard the effects of the free 
molecular static instability. However, if the entry 
spin rate is too large, the gyroscopic stiffness could 
overwhelm the aerodynamic stability in the contin- 
uum regime. This would lead to large angles-of- 
attack during peak heating. After a detailed investi- 
gation performing numerous six-DOF entry analyses 
for a variety of spin rates, an entry spin rate of 13.5 
rpm is selected. This spin rate adequately reduces 
the high altitude angle-of-attack excursions, yet 
avoids any attitude concerns during peak heating. 
Additionally, a 13.5 rpm spin rate affords sufficient 
margin (in the angle-of-attack excursions) to ac- 
commodate off-nominal conditions which may be 
present during the entry (as confirmed by the Monte 
Carlo analysis presented later). 

The transonic/subsonic instability is addressed 
by deploying a supersonic drogue chute prior to 
main parachute deployment. The drogue chute 
serves as a stabilizing mechanism for the SRC until 
main parachute deployment. The drogue chute size 
and deployment Mach number are constrained by 
the need to prevent excessive drift, which could lead 
to a landing footprint beyond the proposed UTTR 
site. However, the drogue chute size must provide 
sufficient area to stabilize the SRC. Furthermore, the 
deployment Mach number must be outside the dy- 
namic instability region near transonic speeds to 
avoid the possibility of large angles-of-attack. From 
spin tunnel tests, the drogue chute is sized to 
0.828 m in diameter to provide ample area for stabi- 
lizing the SRC. 5 Numerous six- and three-DOF 
analyses of the entry reveal that drogue deployment 
at Mach 1 .4 avoids excessive drift concerns (as con- 
firmed by the Monte Carlo results presented later). 


Adoption of these changes into the mission re- 
quired modification of the entire terminal descent 
procedure of the entry. A new deployment algo- 
rithm, consisting of a g-switch (i.e. gravity-switch) 
and two timers, is utilized for deployment of the 
drogue and main parachutes. Previously, only a 
baro-switch was needed for deploying the main 
parachute (diameter = 8.2 m). Figure 5 shows the 
modified nominal entry profile, with the terminal 
descent sequence highlighted. The g-switch is trig- 
gered after sensing 3 g’s, at which point, the drogue 
timer is initiated. After 15.04 sec, the drogue chute 
is deployed, initiating the main timer. After 350.6 
seconds, the main parachute is deployed. This new 
nominal entry sequence is sufficiently robust to ac- 
commodate off-nominal conditions during the entry 
(as confirmed by the Monte Carlo analysis presented 
below). 



Figure 5. Modified nominal Stardust entry sequence. 
Table 3. Nominal Mass Properties of the SRC 


Mass, kg 45.8 

Langley Axis System, 

Center of gravity, m 
Along spin axis (x-direction, 

from nose) 0.283 (X cg /D = 0.3487) 

Off spin axis (y-direction)... -0.00127 
Off spin axis (z-direction) ... -0.00129 


I xx , kg-m 2 (spin axis) 

2.4478 

I VV , kg-m 2 

1.7045 

I zz , kg-m 2 

1.8785 

Ixy, kg-m 2 

0.000816 

Ixz, kg-m 2 

0.000692 

lyz, kg-m 2 

0.000362 
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Trajectory calculations are repeated for the 
modified entry profile using the final, launch con- 
figuration mass properties for the SRC (Table 3). 
The flight characteristics of the modified nominal 
entry are shown in Figs. 6-8. During the entry, the 
SRC aerodynamically decelerates from 12.6 km/s 
(planet-relative velocity) to subsonic speeds. The 
maximum deceleration experienced by the SRC 
during the descent is 33.4 g’s. 



Figure 6. Modified nominal Stardust attitude profile. 



Figure 7. Modified nominal mission entry sequence. 



Figure 8. Modified nominal mission parachute 
deployment sequence. 


Recall that the SRC is statically unstable in the 
free molecular regime. The higher spin rate only 


delays the effect of the static instability so that the 
SRC can traverse through the transitional regime to 
the stable continuum regime without experiencing a 
large increase in the total angle-of-attack. As seen in 
Fig. 6, the total angle-of-attack pitches up to ap- 
proximately 7.0 deg in the transitional regime before 
reducing to less than 2.0 deg near peak heating 
(which occurs around Mach 35.2). Reference 13 
describes the heating environment encountered dur- 
ing the entry. 

As the SRC descends, the static margin de- 
creases near Mach 12 to produce a new trim point. 
Consequently, since the SRC has a non-zero c.g. off- 
set from the spin axis, an increase in a T is observed 
from a mean a T of approximately 1.5 deg near Mach 
12 to approximately a mean a T of 2.0 deg near 
Mach 2. In transitioning to a new trim point, attitude 
rates induce an overshoot in a T (peaking around 
Mach 8) before receding around Mach 2. As the 
SRC approaches transonic speeds, the dynamic in- 
stability drives another increase in a T until drogue 
chute deployment. 

Beginning at Mach 1.4 (approximately 34 km 
altitude), the terminal descent phase of the entry 
begins, which slows the SRC down to approximately 
4 m/s prior to landing. Figure 8 shows the nominal 
altitudes of the drogue and main parachute deploy- 
ments. 

Monte Carlo Dispersion Analysis 

Independent Uncertainty Effects 

Before a combination of off-nominal conditions 
are examined, a sensitivity analysis is first per- 
formed to identify the mission uncertainties which 
have the greatest impact on the overall landing foot- 
print. Each of the 41 mission uncertainties are varied 
independently at their respective ±3-a (maxi- 
mum/minimum) variance. Figure 9 shows the re- 
sulting total downrange obtained from the largest 
contributors to the overall landing footprint. Those 
mission uncertainties which are not depicted lead to 
downrange dispersions less than 0.5 km. 

The mission uncertainties shown in Fig. 9 can 
be grouped into two categories: large contributors 
(mission uncertainties 1-4) and small contributors 
(mission uncertainties 5-9). The first group, con- 
taining initial state vector and atmospheric wind and 
density uncertainties, contribute on the order of 20- 
25 km each to the landing footprint size. Again, 
since the atmospheric winds have a significant im- 
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pact on the downrange due to parachute drift, the 
selection of an appropriate drogue chute size and 
deployment Mach number is critical. The second 
group, containing uncertainties in bus separation 
velocity and aerodynamic drag, produce downrange 
dispersions of approximately 1-5 km each. 


B +3-o uncertainty 
IH1 -3-a uncertainty 

Mission Uncertainty 

1 East-West wind 

2 Density 

3 North-South wind 

4 Initial state vector 

5 Vertical wind 

6 C^: hypersonic-continuum (±4%) 

7 Sep. AV: x-body axis (±0.0508 m/s) 

8 C^: supersonic-continuum (±10%) 

9 Sep. AV: z-body axis (±0.0097 m/s) 


Mission uncertainty 



Figure 9. Significant contributors to the total range 
dispersion. (3-0 variance shown in parenthesis.) 

Multiple Uncertainty Effects 

To determine the robustness of the Stardust 
SRC entry profile, off-nominal conditions are simu- 
lated to address uncertainties which may arise dur- 
ing the descent. The impact of multiple uncertainties 
occurring simultaneously is ascertained by per- 
forming a Monte Carlo dispersion analysis. Over 
3000 random trajectories are simulated to assure 
proper Gaussian or uniform distributions for the 41 
mission uncertainties identified. 



Figure 10. Distribution of total angle-of-attack at 
atmospheric interface resulting from 3000 Monte 
Carlo simulation cases. 



The statistical results from the 3000 Monte 
Carlo simulations are displayed in Figs. 10-17. Fig- 
ures 10-12 show the distribution of the total angle- 
of-attack at three discrete locations during the early 
phase of the mission: at atmospheric interface, in the 
transitional regime, and at peak heating. At atmos- 
pheric interface, the statistical mean total angle-of- 
attack of the 3000 Monte Carlo cases is 2.6 deg. The 
maximum a T observed is around 5.3 deg (which is 
below the mission constraint of 10 deg). In the tran- 
sitional regime, the total angle-of-attack does in- 
crease from atmospheric interface due to the free 
molecular instability. The mean a T is 8.5 deg, and 
the maximum a T observed is 31.4 deg. The higher 
spin rate prevents continued growth in the total an- 
gle-of-attack; so that, by peak heating (where the 
SRC is stable), the mean a T damps to 2.4 deg as 
seen in Fig. 12. The maximum a T observed at peak 
heating is 9.1 deg, which is below the mission con- 
straint of 10 deg. 


Figure 11. Distribution of total angle-of-attack in 
transitional regime resulting from 3000 Monte 
Carlo simulation cases. 



Figure 12. Distribution of total angle-of-attack at 
peak heating resulting from 3000 Monte Carlo 
simulation cases. 
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Figures 13-15 show the distribution of the 
drogue and main parachute deployment conditions. 
The mean Mach number at drogue chute deployment 
is 1.37, as seen in Fig. 13. The minimum deploy- 
ment Mach number encountered is 1.08, which is 
high enough to avoid the significant effects of the 
transonic dynamic instability. The corresponding 
mean total angle-of-attack at drogue chute deploy- 
ment (see Fig. 14) is 11.3 deg. The maximum a T 
observed is 30.5 deg, which slightly exceeds the 
mission constraint of 30 deg. Analysis of this angle 
of attack violation does not pose a concern for a suc- 
cessful drogue deployment. Figure 15 shows the 
distribution of the main parachute deployment alti- 
tude. The mean deployment altitude is 3.08 km, with 
a minimum occurring at 2.26 km. 



Figure 13. Distribution of Mach number at drogue 
chute deployment resulting from 3000 Monte Carlo 
simulation cases. 



a T , deg 


Figure 14. Distribution of total angle-of-attack at 
drogue chute deployment resulting from 3000 Monte 
Carlo simulation cases. 



Altitude, m 


Figure 15. Distribution of altitude at main chute 
deployment resulting from 3000 Monte Carlo 
simulation cases. 


Figures 16 and 17 show the resulting distribu- 
tions in downrange and crossrange at landing for the 
3000 Monte Carlo cases, respectively. The minimum 
downrange is -33.1 km (short) from the nominal 
landing point, whereas the maximum downrange is 
39.2 km (long). The maximum crossrange obtained 
is 13.2 km from the nominal landing point. The re- 
sulting 3-o ellipse has a major axis of 60.8 km 
(-30.6 short, 30.2 long) in downrange and a minor 
axis of 19.9 km in crossrange. This footprint is 
within the UTTR site. Within the assumptions of the 
present analysis, a 99.7 percent probability exists 
that the SRC will land within this 3-a footprint el- 
lipse. Figure 18 shows the landing location of all 
3000 Monte Carlo cases. Table 4 summaries these 
results. 



Figure 16. Downrange distribution at landing 
resulting from 3000 Monte Carlo simulation cases. 
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Figure 1 7. Crossrange distribution at landing 
resulting from 3000 Monte Carlo simulation cases. 
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Figure 18. Landing range dispersion resulting from 
3000 Monte Carlo simulation cases. 

SUMMARY 


Table 4. Summary of Monte Carlo Analysis 



Mean Min. 

Max. 

3-ct 

Attitude dispersion 

Atmospheric interface 




ap, deg 

Transitional regime 

...2.6 

0.5 

5.3 

2.7 

ap, deg 

Peak heating 

....8.5 

3.2 

31.4 

9.8 

ap, deg 

...2.4 

1.1 

9.1 

2.6 

Drogue chute deployment 




ap, deg 

7 5 

0.1 

30.5 

15.5 

Landing dispersion 

Landing 





downrange, km.... 

..-0.6 

-33.1 

39.2 

30.2 (long) 
-30.6 (shorf 

Landing 





crossrange, km 

..-0.1 

-13.2 

11.6 

9.9 

Total range, km 

89 

0.1 

39.2 

17.6 


The solution selected to address the high alti- 
tude instability is to increase the SRC entry spin rate 
to 13.5 rpm. The higher spin rate, although not 
eliminating the instability, increases the gyroscopic 
stiffness of the SRC thereby retarding the effects of 
the free molecular static instability. To address the 
transonic/subsonic instability, a drogue chute (hav- 
ing a diameter of 0.828 m) is added and a deploy- 
ment algorithm based on an gravity-switch activated 
timer resulting in a Mach 1.4 deployment is defined. 
The drogue chute serves to stabilize the SRC until 
main parachute deployment. 


A six-degree of freedom analysis of the nominal 
Stardust SRC entry reveals that two aerodynamic 
instabilities result in unacceptable capsule dynamics 
during the descent using the original entry strategy. 
The first instability resides in the high altitude free 
molecular regime, while the second appears toward 
the end of the entry in the transonic/subsonic flow 
regime. These instabilities, if not eliminated or at 
least suppressed, could lead to mission failure. In the 
free molecular regime, a static instability exists 
which produces large excursions in the angle-of- 
attack (approaching 70 deg) early in the nominal 
entry profile. If off-nominal attitudes or attitude 
rates exist at atmospheric interface, a backward en- 
try is very possible. In the transonic/subsonic re- 
gime, a dynamic instability is present which could 
induce a tumbling motion prior to parachute de- 
ployment. 


For this mission, 41 potential uncertainties were 
identified which could affect the entry. Initial state 
vector and atmospheric property (density, and 
North-South and East-West wind components) un- 
certainties were found to produce the greatest down- 
range dispersions on the order of 20-25 km each. 
Uncertainties from bus separation and aerodynamics 
produced dispersion between 1-5 km each. All other 
uncertainties resulted in dispersion less than 0.5 km. 

A Monte Carlo analysis of over 3000 off- 
nominal trajectories shows that the SRC attitude 
near peak heating and drogue chute deployment to 
be within Stardust mission limits. The resulting 3-o 
landing footprint obtained was 60.8 km (-30.6 short, 
30.2 long) in downrange and 19.9 km in crossrange 
(which is within the Utah Test and Training Range 
boundaries). Within the assumptions of the present 
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study, a 99.7 percent probability exists that the Star- 
dust SRC will land within this 3-a ellipse. 

The instabilities in the Stardust SRC were re- 
vealed too late in the design process to affect the 
design of the capsule. If identified earlier, these 
types of instabilities could be eliminated by consid- 
ering alternative capsule configurations that avoid 
the need for corrective measures late in a program. 
Therefore, a case is made for including six-degree of 
freedom entry trajectory analyses early in the con- 
ceptual design phase. 

Finally, the resolution of the Stardust SRC in- 
stabilities relies heavily on the validity of the Monte 
Carlo analysis. Increased dependence on entry 
simulations for mission success places considerable 
importance on selecting appropriate uncertainties. 
As confidence increases in the analysis accuracy, 
cheaper and/or higher performance entry systems 
can be selected for future missions. As an example, 
for this mission, an inherently unstable capsule can 
still be flown without violating any mission con- 
straints. 
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